CERN-TH/97-246 



UNSTABLE SYSTEMS 
IN 

RELATIVISTIC QUANTUM FIELD THEORY 



L. Maiani\ M. Testa^'^ 
^ Dipartimento di Fisica, Universita di Roma "La Sapienza" 
Sezione INFN di Roma 
Rle A. Moro 2, 00185 Roma, Raly 
2 Theory Division, CERN, 1211 Geneva 23, Switzerland*. 



Abstract 

We show how the state of an unstable particle can be defined in terms of stable 
asymptotic states. This general definition is used to discuss and to solve some old 
problems connected with the short-time and large-time behaviour of the non-decay 
amplitude. 



CERN-TH/97-246 
September 1997 



* Address from Sept 1st, 1997 to August 31st, 1998. 



1 Introduction 



The definition of unstable states in quantum mechanics is notoriously difficult Q , , Q . 
The most usual approach is based on the Breit-Wigner approximation of scattering 
amplitudes @ in which the unstable state is associated with a factorizable simple pole in 
the complex s-channel variable, s being the total momentum squared in s-channel. This 
approach, however, is not suitable to describe the detailed space-time behaviour of the 
unstable system. In particular it underscores the deviations from simple exponential 
behaviour at very short and very large times, which follow from general arguments Q. 

The short-time behaviour of an unstable state weakly coupled to the final decay 
channels has been investigated recently in Ref.[^. The result has been to underline 
an unexpected singularity in the short-time expansion of the non-decay amplitude, 
which makes the deviations from exponential behaviour in this region depend from 
the formation process of the unstable state itself. The short-time singularity and the 
corresponding dependence from the formation condition of the unstable state seems to 
eliminate the so-called Zeno paradox (Q- |ll3[| ) at least for systems described in terms 
of Relativistic Quantum Field Theory. 

In this paper, we present a general, relativistically invariant, definition of the un- 
stable state which allows to elucidate the space-time behaviour of the non-decay ampli- 
tude. The definition follows closely the experimental procedure to measure the lifetime 
of weakly decaying particles (i.e. the so-called impact parameter distribution). 

For weakly unstable particles, where perturbation theory can be used, we find that 
the general quantum-mechanical expression of the decay rate per unit time (the Fermi 
Golden Rule) has to be modified at short times in a well defined way. In addition, we 
recover the results of the previous analysis Q, with a much clearer characterization of 
the non- universal dependence upon the wave-packet shape of the initial particles from 
which the unstable state was formed. 

The space-time behaviour of the non-decay amplitude is specified in terms of S- 
matrix amplitudes. Adopting for the latter the Breit-Wigner form with a constant 
width gives, for times larger than the formation time, the familiar exponential be- 
haviour. 

Finally, on the basis of unitarity arguments, one can estimate the very large time 
behaviour where again deviations from the exponential form are to be expected, in the 
form of a power-law time dependence. It is worth stressing that, in this case, the time 
behaviour is again a non-universal one, related to the form of the initial state wave 
packets. 

Although not strictly related to the central arguments of the present paper, we also 
give, in the last section, a simple and general derivation of the so-called Wigner's delay 



relation|15] for scattering systems, because it fits very naturally within the general 



spatio-temporal description of the scattering process used in this investigation. 

2 Constructing the Unstable State from Scratch 

We define as usual the S- and the T- matrices, according to: 

Sa/3 = {a; out I /3; in) = (a; in\ S \(3; in) = {a; out\ S \(3; out) 



1 



S = I + iT = ^ \6; in) {6; out\ 



^5 ' 

TT+ = T+T = 2ImT (1) 

where a, /3 and 5 denote a suitable set of free particle quantum numbers, e.g. momenta 
and spin components. However, to describe the space-time evolution of the scattering 
process it is necessary to introduce wave-packtes for the initial states. We consider, 
therefore, a two-particle state in the far past: 



\in) 



We consider, for simplicity, equal mass, spinless, particles and we set ourselves in 
the center of mass (c.o.m.) frame of reference. We are interested in a situation in 
which the incoming particles can create a resonant state. We therefore choose the 
wave packets so that f{p) is peaked around some momentum p^^^ and g{p) around 
momentum —p such that the c.o.m. energy is about equal to the resonance mass 



i-res 

M: 



2u;ip) « « M (3) 



Wave packets are choosen so as to represent very distantly localized particles, at 
some large, negative time t = —T, and to overlap around the origin of coordinates, at 
time t = 0. Neglecting long- range forces, which can be dealt with separately, wave 
packets evolve in time freely for all times before collision, up to the last 10~^^ seconds 
or so. We consider the state: 

- iT \in) = (/ - 5) \in) = \in) - \out) = \R) (4) 

The state \out), in Eq.[Q], represents a two particle state in the distant future in 
the wave packets / and g, such that they were overlapping around the origin at t = 0. 
Thus the state vector \R) defined by Eq.Q is the initial state minus the state where 
nothing happens, i.e. it represents the state of the products of the collision which has 
taken place around time t = 0. 

We consider next the amplitude: 

A{t) = {in\ T+e''^^T \in) = {Rt \ R) (5) 

A{t) represents the overlap of the collision state, when the collision has taken place 
at time t = 0, \R), with the collision state when the collision has taken place at a later 
time t (the state {in \ iT~^e~'^^^ = {Rt\). If the collision goes through the formation of 
an unstable state, A{t) is, apart from a normalization factor, the amplitude for this 
state to have remained unchanged during time t, i.e. the non-decay amplitude, so that: 



■^non- decay ) — ^^Q'J 



Pnon— decay (^) 



A{t) 



A{0) 

Eq.p is our basic starting point. We write: 
T 



2 

(6) 



Pi'P2'^^) = (2vr)^(5(^)(p3 -hp4 -Pi -p2)T{s,t) (7) 
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where s and t are the usual Mandelstam variables. Using the unitarity relation, Eq.[^ 
can then be rewritten as: 



X (27r)45W {p3 +Pi-pi- P2)2ImT{s, t) = (8) 
^e-^(^fi+^f2)*(27r)^5W(p3 +P4 - Pi -p2)2ImT{s,t)^ 



where brackets indicate folding with wave packets. Eq.|Q can be rewritten as: 

A{t) = e-'"^* f dxe-'^^'Fix) (9) 



with F{x) a positive definite function with a limited-from-below support (correspond- 
ing to the combined support of ImT and of the wave packets): 

F{x) = ((27r)45(4) (p3 +PA-P1- P2)5{ujq, + cOg, - M - x)2ImT{s, t)) (10) 



We add a few comments. 

• The definition of the non-decay probability, Eqs.|5|,[p, is given in term of S- 
matrix elements, i.e. it involves only the stable asymptotic states out of which the 
unstable state is formed (and in which it will, eventually, decay). No assumption 
is made, in particular, about the existence of a state representing the unstable 
particle itself (at time t = 0), which is indeed a very questionable assumption, 
certainly not valid beyond perturbation theory. If we had assumed that, we could 
have written: 

A{t) = {R\ e-^* \R) = e-^"* |(n I R)? (11) 

The amplitude A{t) in Eq.|5|, as shown by Eqs.|^ and |]l^, shares with the naive 
amplitude, Eq.||Tl|, the property of being the Fourier-transform of a positive def- 
inite function with a limited-from-below support. This property, first stressed by 
L.A. KhalfinQ for the naive amplitude Eq.|ll], gives rise to the non-exponential 



behaviour af A{t) at very large times, as discussed in Section y. 

The above definition of A{t) makes sense only for times t which are larger than 
either the overlap time of the initial wave packets or the characteristic decay time 
of the background processes. In turn, this implies that the definition is useful only 
if the lifetime of the resonance is larger than either these characteristic times. 
Since: 

^overlap ^ ^ ~KE ^ Jd ^^^^ 



the first condition gives: 



r = l«— i «M (13) 

\ overlap 



Similarly, background processes (e.g. box- diagram contributions to the scattering 
process) decay in time with the only time-scale available, that is: 
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so that, in conclusion, our definition is suitable for: 

r = -«M (15) 

T 

as intuitively expected. 

• If we transform from the c.o.m. to another frame of reference, the time translation 
gives rise to a non vanishing space translation. The amplitude A{t) transforms 
into the probability amplitude for the final particles to originate at a non- van- 
ishing distance from the collision point. P{t) corresponds, in this case, to the 
so-called impact parameter distribution of the decay products of an unstable 
particle produced with non- vanishing velocity. 



3 Perturbation Theory 



In this section we show how the previous general scheme works when perturbation 
theory is reliable. In particular we will discuss the derivation of the Fermi Golden Rule 
and its modifications at short times, related to the singularities inherent to Relativistic 
Quantum Field Theory. 

We shall study, for definiteness, a system consisting of an unstable scalar particle 
(p, with mass M, decaying into two scalar particles ipi and ip2 with masses mi and m2 
respectively. We take a decay hamiltonian, Hj, of the simple (super-renormalizable) 
form: 

Hi = g I (fxcl)iJii,2 (16) 



The unstable particle cj) contributes to the elastic scattering amplitude of the two 
particles V'l and ^^2 as: 

Sfi = {p3,Pr,out I pi,p2;in) = dfi + ={2it)^5^'^\pi + P2 - Ps - Pi) (17) 



where 

P2 _ M2 ie ^ (P2 _M'^ + i£) ^ ' (P2 _ M2 



and: 

P =Pl+P2=PZ+P4: (19) 

m(p2) = sm"^ + J dx (0| T{O{x)O{0)) |0) exp{iPx) (20) 

0{x)=gMx)M^) (21) 
The counterterm SM"^ in Eq.[^ is adjusted so that: 

i?en(M2) = (22) 
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and therefore M is the renormahzed mass of the resonant state. The behaviour of 
n(P^) around M is then: 

n(p2) p^-^^.J^P'^ - M'^)ReW [M'^) + ilmn(p2) (23) 

where ReW^M"^) is an ultraviolet divergent parameter related to the wave function 
renormalization of the unstable particle propagator, which is to be removed by renor- 
malization. It will be clear in a moment that this term does not affect the non-decay 
amplitude of the unstable state. As for /mn(P^), it has the expression: 

/mn(p2) = 1 ^ |(o| 0(0) |n)p (2^)45W(P - P„) (24) 

n 

The overlap A{t) defined in Eq.||5|, then becomes: 

A{t) = {Rt I R) = I r{p,)9*{p.)f{pi)9{P2) ^ (25) 

x(27r)^5(^)(pi +P2 -P3 -P4)exp[-i(wpj + ujp^)t]{-iMfi - {iMji)*} 
We have: 



iMf. - {zMf.r = {g?{2^)5{P' - M') - 2g'lm f ^^^J^i^L^ 



and: 



Im 



n(p2) 



(p2 _ M2 + ie)2 



so that: 



(26) 



/ n(p2) \ /mn(p2) 9. . 9 9. 

( (P2_M.;,,). J = (p2_V)2 - 2nReU\M^)5iP^ - M^) (28) 



Eq.||25[] then becomes: 

{Rt I i?) = / nP^)9*{p.)f{pi)9{P2) ^ (29) 
J (21^)° ^2LOp.p.ijjp^2ujp-^2ujp^ 

X (27r)^(5^^) (pi + P2 - P3 - P4) exp[-i(u;p, + Wpjt] x 
X <^ g^{\ + 2i2en'(M2))27r(5(p2 - M^) - 2^^ 



(P2 _ Af 2)2 j 

For simplicity we choose to work in the narrow wave packet approximation. This 
means that we parametrize the momenta as: 



P=P ■ + (30) 

L-t —res % — ^ ' 

and neglect terms quadratic in ki. 



Eq.||29[] can therefore be rewritten as: 

{Rt\ R) = {Rt\ R)-^ + {Rt\ R)^ (31) 

with: 

x{27r)'^6''^\ki + k2-k3-k^)6[2Mv-(k^-k2)]6[v-{k3-k^)] (32) 

and: 



{Rt I R), = -2g^eM-^Mt] f /*(P3)g*(P4)/(pi)g(P2) ^ 3 

^ ' ^2 y PL I J ^27T)(^^2LOp,2LUp,2LJp,2LOp, ^ ' 

x(2vr)V^)(,, -,3 -P4)exp-.[,. (A,, - .,)]tp^^^^ 



Using Eq.[|2^, Eq.[33| becomes: 

{Rt I = 

= -2vr<,2expHMt] ^ K0| 0(0) \n)\' ('^-f^^'^iPn f'^j^^^^^:^^ x (34) 

^ / /*(P3)g*(P4)/(pl)g(P2) ^4.(4)/ , AxrA/T i ^; ^ F 1 

J {2-k)'° y/2ujp.j2u)p^2u)p-,2u]p2 
We can now define q, oq, ai and (3i through: 

{Rt \R)]^ = 2TTg^a exp -iMt = {2TTg^){ao + qi) exp -iMt (35) 

and 

{Rt \R)2 = 27rg^pi{t) exp -iMt (36) 
where a is independent of t, oq is order in g^, while ai and /3i(t) are first order in 

We finahy have for the properly normalized non-decay probability defined in Eq. |p : 

_ \{Rt\R)\' _ \{Rt\R), + {Rt\R)/ ^ 

^non— decay \J') 



\{RQ\R)Y \{RQ\R)^ + {R0\R)2\ 



1 + 2Re[ ao 



1 + 2Re[ — ^-^] (37) 



ai+ft(0) l 

If, for simplicity, we take real momentum- space wave functions, we get: 

|2 



Pnon-decayit) = 1 " ^ E ^^^^^ (2.)35(3) (£ J sin^l^^^t] i/(i?„ - M) 



(38) 



where: 



- M] = (39) 
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is a function which provides a cutoff in energy at small times t, when the energy 
conservation is not yet active. In fact, due to the finite spread in energy of the wave 
packets, H{x) rapidly vanishes for large x, while it goes to unity for vanishing x: 

H[0] = l (40) 

Eq.||3^ reproduces the result first found in Ref.|^]. The " form- factor" , H{En — M), 
provides the cut-off to the sum over intermediate states which is generally needed to 
cope with the singular behavior at very small times, as further discussed in the next 
section. At larger times, t, on the other hand: 

En-M. 



,sin[^^^t] 



^ 2TTt6{En - M) (41) 

and Eq.||3^, in virtue of Eq.||40[|, reduces to the usual Golden Rule formula for the 
non-decay probability of an unstable particle: 

Pnon-decay{t) = 1 " ^ E ^(0) \n)\^ (27r)^ ^^^^^ (£„)5[^n " M] (42) 



n 



4 The Nature of Short-Time Singularities 

As discussed in Ref.|^, the use of perturbation theory without taking into account 
the formation time of the resonant state, would lead to an expression for Pnon- decay {t) 
similar to the one given by Eq.[^ with H{En — M) = 1, i.e.: 

P„.-.„.(*) ^ 1 - A X: K«^p,,3„3,,£^) (43) 

It was shown in Ref.[^ that the expansion of Eq.|43| in powers of t is usually marred 
by meaningless ultraviolet divergencies. In this section we discuss in more detail the 
nature of these short-time singularities. This discussion clarifies the nature of the 
assumptions needed to derive the so called Zeno paradox (0-[|l^]). In fact the Zeno 
paradox strongly depends on the quadratic short-time behaviour, usually inferred from 
the (highly formal) argument: 

Pnon— decay {i) l^-^l^ ' l-^^l 

= 1 - t'^{< P \ H"^ \ P >- \ < P \ H \ P > \ ^) + ... = 

= 1 - AE'^t'^ + ... (44) 



Expanding the explicit expression of Pnon- decay [t) given by Eq.[Q in powers of 
we have: 



Pnon-decay{t) = 1 " ^ E ^(0) \n)\\27,)H^^HP^) = 

= 1-^/ dx(0|O(x,0)O(0)|0) (45) 
Eq.E^ shows that in this way we get, indeed, formally, a short-time behaviour. 



as in Eq.[44|. The problem is that the operator product appearing in Eq.|45] is not 
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integrable. In fact, apart from possible logs, we have from the Operator Product 
Expansion: 

0(x, 0)0(0)^^^^/ + ... (46) 

where do denotes the dimension of the operator O. In the present example c?o = 2 
and Eq.[^5| suffers from a linear ultraviolet divergence. The situation worsens for more 
singular (higher dimensional) decay hamiltonians, as discussed in Ref.Q. 

These considerations clearly show that there is a short-time ultraviolet singularity 
which: 

• makes Ai?^ divergent in Eq. 



• is smeared, in the particle survival amplitude Eq.|38|, over time scales provided by 
the formation mechanism of the unstable state (the overlap time of the incident 
wave-packets) . 

It should be clear that no form factors (in the case of decays involving hadrons, as 
e.g. proton decay) can cure this divergence which is quite similar to the one measured 
in deep-inelastic- scattering, due to the singular product of two currents. 

In Eq.[38|, we can interpret the presence of H{En — M) as a cut-off on the energy 



of detected final decay products (a reasonable requirement for a measuring apparatus) , 
which, by the optical theorem, produces an effective smearing on the particle survival 
amplitude. 

We stress again the conclusions of Ref . Q concerning the fallacy of the finiteness of 
AE'^ in Eq.[^], on which many of the papers dealing with the Zeno paradox are based. 



5 The Exponential Decay 



Higher orders in perturbation theory give rise to higher powers of the time, t, in Eq.|42|. 
For times of the order of the lifetime r we need to sum up at least the leading terms, 
to get a meaningful result. 

The linear behaviour in time, in Eq.[^], arises from the singular behaviour in 
(s — M^) of the self-energy insertion. Thus, higher powers in t shall correspond to 
repeated self-energy insertions in the lowest-order particle propagator. 

More precisely, the leading singular behaviour is given by insertion of the imaginary 
part of the self-energy evaluated at s = M^. Inserting the successive terms of the 
expansion of /mil around s = would give rise to less singular terms in (s — M^), 
i.e. to a subleading behaviour in time. The same applies to the insertion of -Rell, since 
it vanishes at s = by mass renormalization. 

Using the above considerations, we find that the sum of the leading terms in time 
is given by the Breit-Wigner propagator: 

Ait) = (e-*('^-i+-''2)*(27r)45(^)(p3 +P4 -Pi -P2)2/m[r(s,i)]w) (47) 



= .-M2 + ™ ^^^^ 
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A{t) = e"^*" / dxe-^^* ^2, Gix) (49) 



so that (compare with Eq.||9[]): 

2MT 

4M2(x2 + ^) 

with G{x) defined by: 

G{x) = ((27r)45W (p3 +P4-P1- P2)6{iOp, +uJp,-M-x)) (50) 
In the narrow wave packet approximation, Eq.||3^, we have: 

ujp^+ojp^^ M + y_- (ti - ^2) (51) 
The region where G{x) ~ constant ~ G(0) is thus hmited by: 

V 1 

Xcrit ~ vAp « — PS — >> r (52) 

lAX overlap 



Since the integral in Eq.|49| is dominated by the region x ~ F, we may take G{x) 
as a constant, to a good approximation. Extending, furthermore, the integration range 
to ±00, we find: 

A(t) = I^e-"'.-¥ (53) 
that is, in conclusion, a pure exponential behaviour for the non-decay probability. 

6 Very Large Time Behaviour 

The exponential behaviour displayed in Eq.[^] is the result of an approximation 
In fact we kow from the Riemann-Lebesgue lemma that the asymptotic behaviour of a 
Fourier transform integral, as the one appearing in Eq.j^, is determined by the points 
where F(x) or some of its derivatives are singular. If F{x) where continuous together 
with all its derivatives, the non-decay amplitude Anon-decay{t) whould vanish faster 
than any power of i. The presence of any singularity makes Anon- decay (t) vanish not 
faster than some power of j. From the expression of F{x), Eq.||l^, it is clear that 
singularities are necessarily present which originate from two possible sources: 

1. singularities due to the unitarity of the 5'-matrix; 

2. singularities of the wave function of the initial state. 

As for 1), unitarity requires that a resonance pole be located on the second sheet of 
the complex energy plane, thus implying the existence of at least one branch singularity 
of the amplitude T{s,t) in Eq.pO|]. 

The presence of this singularity will, however, be ineffective for the large t behaviour 
of Anon-decay{t), unless the resonancc location is very close to threshold. In fact, as 
can be seen from Eq.[10|, its contribution is, in general, depressed by the narrowness 



of the wave packet of the initial state. 

On the second point very little can be said. The possible singularities due to the 
initial wave function have to do with the details of the experimental preparation of the 
resonant state and are certainly not even under direct control of the experimentalist. 
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Whatever their origin, the presence of these singularities in F{x) make Anon-decay{t) 
behave asymptotically in a power-like way, rather than exponentially. It must be 
remarked, however, that nothing general can be said on the onset time of this power- 
like behaviour because it depends on the detailed structure of the wave function of the 
initial state. 

To get a crude idea of the times where the power law behaviour takes over, we may 
replace Eq.|48| with a form in which the two particle cut is considered: 

-1 



(54) 



s - M2 + iMn(s) 
where, to account for threshold behavour: 

n(.) = e(.-W)[^]2'+ir (55) 

and: 

Ak{sf = s- Ani^ (56) 

but otherwise we neglect any energy dependence of the matrix element. With these 
positions, Eq.|49| becomes: 

A{t) = exp{-2imt) f dxexp{-ixt)[2MU{s)][{s - M^f + M'^Ii{sf]-^G{x) (57) 
Jo 

with 

y^s) = x + 2m (58) 

Neglecting further any variation of the function G{x), which describes the energy 
spread of the initial wave packets, the asymptotic behaviour of A[t)^ as t — > cxd, is 
easily estimated to be: 

A{t)asyn.pt = C[^)ri{Ttr'~'i (59) 

with: 

^ = (t^)(^)'^W / duexp{-iu)u'+'2 (60) 
128 rw^ Jo 

Comparing with the exponential law, Eq.|5^, we see that the power law gets in for a 
critical time , tcrit, which is very large indeed. Numerically: 

Tt^rit ~ 2(1 + -){25 + ln{\ ^^^^ If JJ^ 1}} (61) 

This result amply justifies the fact that no deviations from the exponential be- 
haviour have been observed, until now, for long-lived systems [ 14 1. 



7 Wigner's Delay 

Within the formalism described in the present paper, we can readily recover, in a 
very general way the result, due to Wigner||l5| and usually derived within potential 
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scattering theory, which says that the derivative of scattering phase shift 6i (E) is related 
to the delay time induced by the interaction through: 

d6i{E) 

T = 2-^ (62) 

In order to prove Eq.[^ let us consider the center of mass scattering of two spinless 
particles in a state of given orbital angular momentum Z, under the inelastic threshold. 
The incoming state (a wave packet normalized in a finite volume V) is: 



I/; in) = y ^ j dE f{E) \E, I, m; in) (63) 

where \E , I , m; in) denotes the incoming state of two particles of given total energy E 
and given angular momentum I, m, normalized as: 

(E', I', m'; in \ E, I, m; in) = j^-^5ii>6,mm'^iE - E') (64) 



and 



dE\f{E)\^ = l (65) 



The amplitude to find, after the interaction, the system in a state which corresponds 
to the free propagation of the initial state is: 

A = {f;out\ f;in) = 
(27r)3 



V 



J dE'dEf*{E')f{E){E',l',m';out \ E,l,m;in) 



= J dE\f{E)\^e^p2i5i{E) (66) 

If the system were non interacting, Si{E) = 0, and we would have A = 1. 
In the case of a narrow wave packet, in which f{E) is strongly peaked around a 
given energy Eq, we have: 

A ^ exp[2i6i{Eo) - 2i5[{E^)E^] j dE \f{E)\^ eM'^i5[{E^)E] (67) 

We can now ask what is the probability to find the system, after the interaction, 
in a state which is the free propagation of the initial state delayed by a time T. Such 
a time- translated state is given by: 



\f,T;out)=expiHT\j^^^ J dE f{E)\E,l,m;out) = 

= J dE f{E) exp[iSr] \E, /, m; out) (68) 

where, as usual, H denotes the full hamiltonian of the interacting system. 
We have, in this case: 

A(T) = {f,T-out\f;in)^ 

^ exp[2i5;(So) - 2i<5K^o)^o] / dE|/(£;)|'expiS[25K^o) - T] (69) 
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From Eq.^^ it follows that for T = 26'i{Eo), the probability becomes 1: 

\AiT = 26'i{Eo))\^ = I (70) 



which proves the Wigner delay relation, Eq.|62|. 
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